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Indian Standard 
USE OF PROBABILITY PAPERS 



0. FOREWORD 



0.1 This Indian Standard was adopted by the 
Bureau of Indian Standards on 23 February 1988, 
after the draft finalized by the Qjiiality Control 
and Industrial Statistics Sectional Committee had 
been approved by the Executive Committee. 

0,2 Some commonly used statistical models have 
been discussed in IS : 9300 ( Part 1 )-1979* and 
IS : 9300 ( Part2 )-1979t. However, the selec- 
tion of a suitable model which represents the 
given data is also of major concern. For this 
purpose, two different approaches may be consi^ 
dered, namely, probability plotting method and 
statistical tests. Probability plotting is a subjective 
method as the determination of whether or not 
the data contradicts the assumed model is based 
on visual examination rather than statistical 
calculation. However, the probability plotting 
method is very simple and may provide a great 
deal of useful information in addition to evalua- 
tion of appropriateness of the chosen model. 

0.3 The technique of probability plotting provides 
a pictorial representation of data as well as the 
following: 

a) Evaluation of reasonableness of the assum- 
ed statistical model. 



♦Statistical models for industrial applications: Part 1 
Discrete models. 

fStatistical models for industrial applications: Part 2 
Continuous models. 



1. SCOPE 

1.1 This standard covers the probability papers 
for normal, log-normal, exponential, Weibull, 
and binomial distributions. The procedure for 
estimating parameter (s) of the distribution has 
also been discussed. Each probability paper 
along with estimation of the parameter (s) has 
been illustrated with examples. Limitations on 
the use of probability papers for parameter (s) 
estimation has also been discussed. 

2. TERMINOLOGY 

2.1 Normal Distribution — The probability 
distribution of a continuous random variable X, 
the probability density function of which is: 



b) Estimation of percentiles of the distribu- 
tion, and 

c) Estimation of distribution parameter (s). 

0.4 In this technique of probability plotting, 
special probability graph papers have been desi- 
gned for the distributions. The data is plotted on 
probability graph paper of an assumed distribu- 
tion. The operational procedure for plotting has 
been given in 4 of this standard. If the assumed 
model is correct, plotted points on the graph 
paper will tend to fall in a straight line. If the 
assumed model is not appropriate, the plot will 
not be linear; and the extent and type of disper- 
sion may be examined. Some guidelines for judg- 
ing the extent of departure have also been 
discussed. 

0.5 Once the chosen model appears to be appro- 
priate for the given data, the parameter(s) of the 
model is ( are ) estimated from the probability 
plot. The procedure for estimation of parameter(s) 
has been illustrated with examples. 

0.6 In reporting the result of a test or analysis, 
if the final value, observed or calculated, is to be 
rounded off, it shall be done in accordance with 
IS : 2-1960*. 



•Rules for rounding off numerical values ( revised ). 



/w 



I r - ( ^ - (X )n 



^ < -v < + 00 



2^2 Log-Normal Distribution — The probabi- 
lity distribution of a continuous random variable 
X, the probability density function of which is: 

ff ^ ^ r -(log x-^Y ^ 

< ^ < 00 

2.3 Exponential Distribution — The probabi- 
lity distribution of a continuous random variable 
X, the probability density function of which is: 

f {x) ^ \ exp ( — A a: ), 

X y 0^ and A > 

2.4 Weibull Distribution — The probability 
distribution of a continuous random variable X, 
the probability density function of which is: 

^(*^ = t(x)''' ^^p [ - ( ^/* )'' j. * > 

a and b are positive constants and determine 
shape and scale of the distribution respectively. 



1 
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2.5 Binomial Distribution — The probability 
distribution of a discrete random variable X, such 
that: 

a: = 0, 1, 2, n 

< p <\ 



where 



(:)- 



a: ! { n — X ) I 



3. LIMITATIONS ON USE OF 
PROBABILITY PAPERS FOR 
PARAMETER ESTIMATION 

3.1 The use of probability papers for estimation 
of parameter is essentially based on the estima- 
tion of population quantiles ( or percentiles ) by 
the corresponding sample quantiles. For estimat- 
ing the two parameters in either the normal or 
the log-normal or the WeibuU distribution, one 
requires to estimate two quantiles using the pro- 
bability plots. The choice of these two quantiles 
is not unique unless determined through the 
complicated approach to minimize the generalized 
variance of the estimators over the possible 
choices. Thus one may use different pairs of 
quantiles to get different pairs of estimates from 
the same data for the two parameters in any of 
these three models. Further, the estimators obtain- 
ed by this graphical method are not quite 
efficient compared to the maximum likelihood 
estimators. These two points should be kept in 
mind while using the technique of probability 
plotting. The principal merit of this technique is 
its simplicity, 

4. OPERATIONAL PROCEDURE FOR A 
PROBABILITY PLOT FROM THE DATA 

4.1 The various steps in preparing a probability 
plot for a continuous distribution from a given 
data are as follows. 

4.1.1 Obtain a probability paper, designed for 
the distribution under examination. 



4.1.4 For censored data, n represents the total 
number of samples, irrespective of whether or not 
all their values are known. The axis of the graph 
paper on which the^j's are plotted will be referr- 
ed to as the observational scale, and the axis for 
(2/ — I ) 50/n as the cumulative scale. Whether 
the observational or cumulative scale appears on 
the abscissa differs for ihe various probability 
papers. However, the cumulative scale is easily 
recognized by the fact that its values range from 
a small positive number, such as O'Ol percent 
( or 000 1 ) to some value close to 100 percent 
( or 1-0 ), such as 99*99 percent ( or 0*999 9 ). 

4.1.5 If a straight line appears to fit the data, 
draw such a line on the graph paper. 

4.1.6 If the chosen model is correct, the points 
should cluster around the line although there will 
be some deviations because of random-sampling 
fluctuations. If the plot deviates appreciably from 
linearity, the chosen model does not appropria- 
tely describe the given data. Systematic depar- 
tures are indications that the model is not 
appropriate. The determination of what may be 
considered a straight line is a subjective matter, 
and two people looking at the same plot may 
arrive at different conclusions. However, the 
larger the number of samples and the greater the 
divergence from the assumed model, the easier it 
will be to detect a true departure. 

5. GUIDELINES FOR JUDGING THE 
APPROPRIATENESS OF A PLOT 

5.1 The following are the general guidelines given 
for judging the appropriateness of the assumed 
model and systematic departures from a stiaight 
line. 

5,1.1 The dispersion of the plotted points in 
the tails ( extreme low or extreme high ) will be 
larger than that of the points at the centre of the 
distribution. Thus, the relative linearity of the 
plot near the tails of the distribution will often 
seem poorer than at the centre of the distribution 
even if the correct model has been chosen. 



4.1.2 Rank the observations in non-decreasing 5.1.2 The plotted points are ordered and hence 

order. Thus, if x (1), x (2), , x {n) arc are not independent. Thus, one should not expect 

the original unordered observations, then the them to be randomly scattered about a line. For 

ordered observations, to be denoted as xi, xj, example, the points immediately following a 

, x^y are such that ^i < ^2 < , point above the line are also likely to be above 

< Xfl, the line. Hence, even if the chosen model is 

- - - „, ^ , , i_ i_-i- correct, the plot may consist of a series of sue- 

4.1.3 Plot the ;ri s on the probability paper ^^^^j^^ j^^^^^ ^^^^^ ^^ ^^^3 ^ ^^^^^ ^^^ ^^^^^ 

versus ( 2j— 1 ) 50/n, or ( 2? — I )/2n dependmg the linef 
on whether the marked axis on the probability 

paper refers to the percent or the proportion of 5.1.3 A model can never be proven to be 

observations. It is henceforth, assumed that the appropriate on the basis of sample data. Thus the 

probability paper requires the plot of percent probability plot of a small sample taken from a 

observations. Thus the lowest observation is near-normal distribution will frequently not differ 

plotted versus 50/n; the second lowest versus appreciably from that of a sample from a normal 

150/«, and so on. distribution. 
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5.1.4 An examination of a plot when there is 
a systematic departure from a straight line can 
give some information about the shape of a more 
appropriate model. For example, if the plotted 
lower points are too large and the upper points 
are too small, compared with those expected from 
the chosen model, it indicates a distribution which 
is less peaked than a normal distribution. 

5.1.5 If the values in both tails are too large, 
the data suggests that a model skewed to the right 
should be chosen. If the values in both tails of 
the normal probability plot are too small, a distri- 
bution skewed to the left would be suggested. 

6. NORMAL PROBABILITY PAPER 

6.1 The normal probability paper is given in 
Fig. I. The procedure for estimation of para- 
meters of normal distribution is given in 6.2. 

6.2 Estimation of Parameters — For normal 
distribution, mean and median are equal. There- 
fore, from the probability plot, the 50th percentile 
or median is estimated as follows: 

a) Locate 50 percent on the cumulative distri- 
bution scale, 

b) Find the intersection of the line obtained 
from the plot and the 50 percent line, 
and 

c) Read the corresponding value on the 
observational scale. 

6,2.1 In the similar way, the estimate of popu- 
lation standard deviation ( o" ) is obtained by 
using the fact that for any normal distribution, 
the standard deviation is approximately equal to 
two-fifth of the difference between the 90th and 
the IGth percentiles. 

6.3 Example — Table 1 gives 40 test results of 
initial setting time ( in minutes ) of ordinary port- 
land cement in non-decreasing order. Test for 
the normality and obtain the parameters of 
normal distribution. 

6.3.1 For each i, the value of /»i = ( 2i — 1 ) 
501 n is calculated ( see Table 2 ). The pairs of 
points ( pi, xi ) are plotted on the normal probabi- 
lity paper ( .^ee Fig. I ). It may be seen from this 
plot that the data seem to follow normal distribu- 
tion. 

6.3.2 For obtaining the parameters, the values 
of Xi corresponding to 50th percentile (see Fig. 1 ) 
is 5\*5, which is the mean of the normal distri- 
bution. For estimation of standard deviation, the 
90th and 10th percentiles are corresponding to the 
values 62 5 and 40*5 on the observational scale. 

A 

Thus, a -= 2 ( 62-5 — 40-5 )/5 = 88 minutes. 



TABLE 1 INITIAL SETTING TIME ( MINUTES ) 

( Clause 6,3 ) 
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45 


50 


54 


59 


36 


46 


50 


54 


60 


38 


47 


51 


55 


61 


40 


47 


51 


55 


63 


41 


48 


52 


56 


65 


42 


49 


52 


57 


67 


43 


49 


53 


57 


69 


43 


50 


53 


58 


70 







TABLE 2 


VALUES OF Xi'i 


3 AND Pi's 










( Clause 6.3.1 ) 
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Xi 


Pl 


i 


^1 


Pi 


1 


Xi 


Pi 


(I) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


35 


1-25 


15 


49 


36-25 


29 


56 


71*25 


2 


36 


3*75 


16 


50 


38-75 


30 


57 


73-75 


3 


38 


6*25 


17 


50 


41-25 


31 


57 


76*25 


4 


40 


8*75 


18 


50 


43-75 


32 


58 


78-75 


5 


41 


11-25 


19 


51 


46*25 


33 


59 


81-25 


6 


42 


13-75 


20 


51 


48-75 


34 


60 


83-75 


7 


43 


16-25 


21 


52 


51-25 


35 


61 


86-25 


8 


43 


18-75 


22 


52 


53*75 


36 


63 


88-75 


9 


45 


2r25 


23 


53 


56*25 


37 


65 


91-25 


10 


46 


23-75 


24 


53 


58-75 


38 


67 


93-75 


11 


47 


26-25 


25 


54 


61-25 


39 


69 


96.25 


12 


47 


28-75 


26 


54 


63-75 


40 


70 


98-75 


13 


48 


31-25 


27 


55 


66-25 








14 


49 


33-75 


28 


55 


68-75 









7. LOG-NORMAL PROBABILITY PAP£R 

7.1 The probability paper for the log-normal 
distribution differs from that for the normal distri- 
bution in the sense that the axis for the observa- 
tional scale is logarithmic instead of arithmetic 
( see Fig. 2 ). A log-normal probability plot can 
also be prepared on normal probability paper by 
using the logarithms of the observations and 
plotting these against ( 2? — 1 ) 50/«. 

7.2 Estimation of Parameters — The proce- 
dure for estimating the parameters is similar to 
that for the normal distribution, except that the 
logarithms of the observed values are used. Thus 
for log-normal distribution, mean ( /Lt ) is estimat- 
ed as the logarithm of the plotted 50th percentile, 
and standard deviation ( a ) is estimated as two- 
fifths of the difference between the logarithms of 
the plotted 90th and 10th percentiles. The log- 
normal probability paper is illustrated with the 
help of example given in 7.3. 
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7.3 Example — A new control device was tested 
on 96 diesel locomotives. Whenever a device 
failed, the distance in kilometres was recorded 
and the device was returned to the factory for 
failure analysis. After 1 35 000 km of service each 
unfailed device was removed from its locomotive. 
It was felt that the device could be warranted 
for ICO 000 km. The underlying failure mecha- 
nism suggested a log-normal distribution for time 
to failure. The ordered times for the 37 devices 
that failed are given in Table 3. Test whether the 
data follows log-normal distribution. Obtain its 
parameters and also the probability that the 
device will not fail in the warranted period. 

7.3.1 The data is censored as the values for 
the remaining 59 devices are known only to be 
above 135 000 km. But the value of w, the sample 
size, in the probability plot has to be taken as 
96 ( see 4.1.4 ). The values of x^% with corres- 
ponding probabilities ( p\ ) are given in Table 4. 

7.3.2 The pair of points ( px^ xi ) are plotted 
on the log-normal probability paper ( see Fig. 2 ) 
which suggests that the log-normaJ distribution is 
a reasonable choice since the plotted points fall 
close to the line drawn. 



TABLE 3 


DISTANCE IN MILLION METRES 




TO FAILURE 






( Clause 7.3 ) 




22-5 


69-5 




93-5 


120*0 


37*5 


76-5 




102*5 


122*5 


46-0 


77-0 




107*0 


123*0 


48-5 


78-5 




108*5 


127*5 


51-5 


80*0 




112*5 


131*0 


53'0 


81-5 




113-5 


132-5 


54-5 


82-0 




116*0 


134*0 


57-5 


83-0 




117*0 




66-5 


84*0 




118*5 




68-0 


91-5 




119*0 




7.3.3 From Fig. 2, 


the 


50th 


percentile is for 



;ti = 1 65 000. Therefore, estimate of mean { ^ ) 
= log ;ci = log ( 165 000 )= 12-01. Similarly 
the values of ^i for 90th and 10th percentiles are 
400 000 and 66 000, respectively. Therefore, 

(T = 2 { lege ( 400 000 ) - ,loge ( 66 000 ) }/5 
= 0-72. 



TABLE 4 VALUES OF xi's AND pi's 
(C/a«i« 7.3.1 ) 



% 


XI 


Pl = 


( 2i - 1 ) 50 
96 




x\ p\ 


( 21 - 1 ) 50 
96 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


I 


22*5 




0*5 


20 


91*5 


20*3 


2 


37*5 




1*6 


21 


93.5 


21-3 


3 


46*0 




2*6 


22 


102*5 


22*4 


4 


48*5 




3*6 


23 


107-0 


23*4 


5 


51-5 




4*7 


24 


108*5 


24*5 


6 


53-0 




5'7 


/25 


112*5 


25*6 


7 


54*5 




6*8 


26 


113*5 


26*6 


8 


57-5 




7-8 


27 


116*0 


27*6 


9 


66*5 




8*9 


28 


117*0 


28*6 


10 


68'0 




9-9 


29 


118*5 


29-7 


11 


69*5 




10*9 


30 


U9-0 


30-7 


12 


76-5 




12*0 


31 


120*0 


31*8 


13 


77*0 




13*0 


32 


122*5 


32*8 


14 


78'5 




14*1 


33 


123*0 


33*9 


15 


80*0 




15-1 


34 


127*5 


34*9 


16 


81*5 




16*1 •' 


35 


131-0 


36*0 


17 


82*0 




17-2 


36 


132*5 


37*0 


18 


83*0 




18*2 


37 


134*0 


38*0 


19 


84*0 




19*3 
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7.3.4 From Fig. 2, the value on cumulative scale 
corresponding to ^i = ICO 000 km is 24 percent 
which gives the probability of failure of device in 
less than 100 000 km. Thus, the estimated reli- 
ability ( probability of no failure ) in the first 
100 000 km is 76 percent, 

8. PROBABILITY PAPER FOR 
EXPONENTIAL DISTRIBUTION 

8.1 The probability paper for exponential distri- 
bution is shov^n in Fig 3. This paper is further 
illustrated with the example given in 8.3. 

8;2 Estimation of Parameter — The para- 
meter of exponential distribution ( A ) is estimated 
by the following relation: 



A = 



Zt - ^1 
2 (J. - 71 ) 



where J 1 and jg are two values on the observa- 
tional scale and Zi and -Ca are the corresponding 
values from the line at the top of the probability 
paper. 

8.3 Example — In a chemical analysis, a given 
contaminant can be measured analytically for 
amounts of 01 percent or more. A sample is 
chosen from each of the 16 lots and the percent 
contaminant is found to be more than 1 
percent for 12 of the 16 analysis. Thus, it is a 
censored sample with the four smallest observa- 
tions missing. The data after ordering are shown 
in Table 5. Determine whether the data fits to 
the exponential distribution and estimate its 
parameter. 

8.3.1 The values ofj&j = ( 2i — 1 ) 50 /« have 
been calculated and given in Table 5. Since this 
is a censored sample with the four smallest obser- 
vations missing, the first point of ^i corresponds 
to the fifth ordered sample and is plotted versus 
(2x5 — 1 ) 50/n, the next value is plotted against 
(2x6 - 1 ) 50/n. 

8.3.2 The closeness of the plotted points to a 
straight line indicates that the model appears to 



represent the data reasonably well. The values of 
the censored observations may be estimated by 
determining values on plotted line corresponding 
to ( 2z — 1 ) 50/n for i = 1, 2, 3 and 4, that is, 
the values on the observational scale correspond- 
ing to the values 3*1, 9 4, 15 6 and 21-9 on the 
cumulative scale. The resulting estimates are 
0-000 5, 002 8, 0004 5 and 0*006 5 percent. 

8.3.3 As mentioned in 8.2, the parameter of 
exponential distribution is estimated as: 

A (3-60 — 2-16) 

A 



36-0 



2 ( 005—003) 
9. WEIBULL PROBABILITY PAPER 

9.1 The use oi Weibull probability paper ( see 
Fig. 4 ) and estimation of its parameters are 
explained with the help of an example given 
in 9.2. 

9.2 Example — The number of cycles to failure 
in a life test of 16 carbide inserts are given in 
Table 6. Test if the data follows Weibull distri- 
bution and obtain its parameters. 



TABLE 6 FAILURE TIME FOR CARBIDE INSERTS 


( Clause 9.2 ) 




210 619 


904 


1 256 


271 698 


1 008 


1 417 


396 811 


1 130 


1 500 


402 900 


1 230 


1 621 



9.2.1 The values oi py% have been calculated in 
Table 7. These values of /»j's when drawn against 
^j's on Weibull probability paper ( see Fig. 4 ) 
indicate a straight line, thereby implying that 
Weibull distribution seems to be reasonable repre- 
sentation of the data. 

9.2.2 To obtain the parameters of Weibull 
distribution, first the slope and intercept of the 
line arc obtained as follows; 



Slope . 



Zi 



TABLE 5 PERCENT CONTAMINANT IN CHEMICAL COMPOUND 

( Clause 8.3 ) 
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*1 
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P\ 


Xi 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 to 4 
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less than 
0-010 


11 


^y^ 


0*030 5 


5 


28*1 


0*012 


12 


71*9 


0*036 3 


6 


34-4 


0*012 8 


13 


78*1 


0*040 2 


7 


40*6 


0*013 9 


U 


84-4 


0*050 3 


8 


46-9 


0-017 1 


15 


90-6 


0*068 9 


9 


53*1 


0*023 I 


16 


96-9 


0*096 8 


10 


59-4 


0*025 
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where W^ and Wi are two values to be read 
from the right ordinate, and J^i ^^^ ^i ^r^ ^^^ 
corresponding values ( to W^ and Wi ) to be read 
from the top abscissa. 

Intercept = c = the value of W to be read 
from the right ordinate where the value of >^ ( to 
be read from top abscissa ) is zero. 

The parameters of Weibull distribution, 
namely, a ( shape parameter ) and b ( scale 
parameter ) are then obtained by the following 
relation: 



m 



b = cxp ( —elm ) 





TABLE 7 


VALUES OF XI 


»8 AND ^i»9 










{Clausi9.2A ) 
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XI 


Pl- 


= { 2i - I ) 
X 50/16 
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^1 


pi 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


1 


210 




3'1 


9 


904 


53-1 


2 


271 




9-4 


10 


1 008 


59*4 


3 


396 




15-6 


11 


1 130 


65*6 


4 


402 




21*9 


12 


1 230 


71-9 


5 


619 




28M 


13 


1 256 


78*1 


6 


698 




34-4 


14 


1417 


84'4 


7 


811 




40*6 


15 


1500 


90-6 


8 


900 




H^ • 


16 


1 621 


96-9 



9.2.3 In this example, for W^ = 20 and 
IVj = -~8'0, the corresponding values of <^s and 
Zi are 798 and 2^88. Therefore,^ 



Slope 



m 



- 2 - ( -8-0 ) 

- 798 — 2-88 

c«- 13-5. 



1-96, and 



Intercept 

Therefore, a = 1-96 and 
^=.exp{-^(-13'5)/ 1-96} = 980. 

10. BINOMIAL PROBABILITY PAPER 

10.1 For the construction of binomial probability 
paper, the horizontal and vertical axis are scaled 
so that from origin the point 1 cm is scaled by l^, 

2 cm by 2^ ,dfcmbyflf«, and so on. The 

points on the plane are denoted by { x^y ), where 
X is the scale of the horizontal axis ( abscissa ) , 
and jf that of the vertical axis, x and j denote the 
number of non-defectives and defectives, respecti- 
vely in a sample of size n(x4-j' = n). A quar- 
ter circle with centre at origin and radius of 100 
[ X =^ y ~ 100 ) is drawn ^ see Fig. 5 ) . 

10.2 Procedure for Estituation of Percent 
Defective and GonEdence Interval — Take a 
sample from an unknown population and count 
the number of non-defectives {x) and that of 
defectives (y). Plot the point { x, y ) on the 
binomial paper ( see Fig. 6 ). Connect the point 




Fio. 5 QjaARTER Circle 

{ x,y ) with the origin. Extrapolate it, and read 
the scale of vertical axis at the intersection of the 
line and the quarter circle. This is the point 
estimate of percent defective. Draw two lines 
around the point { x^y ) with a radius of one 
centimetre. Draw two lines which pass through 
the origin and touch the circle at the upper and 
lower circumferences. Read the scale of vertical 
axis for the intersection of two lines and the 
quarter circle. These two values give the confi- 
dence limits for percent defectives. 




Fig. 6 Binomial Probability Paper 

10.3 Example — When the sample size is 30 
and the number of defectives in the sample is 9, 
the point estimate and the confidence limits of 
fraction defective are obtained as 30 percent and 
15-47 percent, respectively ( see Fig. 6 ). 
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